Solitons in isotropic antiferromagnets: beyond a sigma model. 
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Isotropic antiferromagnets shows a reach variety of magnetic solitons with non-trivial static and 
dynamic properties. One-dimensional soliton elementary excitations have a periodic dispersion 
law. For two-dimensional case, planar antiferromagnetic vortices having non-singular macroscopic 
core with the saturated magnetic moment are present. The dynamic properties of these planar 
antiferromagnetic vortex are characterized by presence of a gyroforce 
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I. INTRODUCTION 

Magnetically ordered materials (magnets) are known 
as essentially nonlinear systems and shows a large va- 
riety of localized nonlinear excitations with finite en- 
ergy, or solitons, see Refs. [HIIIIIH It is sufficient to 
note kink-type solitons (domain walls) which destroy 
long range order in one- dimensional systems; magnetic 
vortices, which cause a Berezinskii-Kosterletz-Thouless 
transition in two-dimensional magnets with continu- 
ous degeneration;^ and also two-dimensional localized 
solitons like Belavin-Polyakov solitons p see for review 
Refs. 0. All these solitons were firstly introduced in 
physics of magnets, and the development of soliton con- 
cept for this particular region of physics is believed to be 
important for modern nonlinear general physics of con- 
densed matter as well as for field models of high-energy 
physics, see Ref. H. 

To date, solitons in Heisenberg ferromagnets, whose 
dynamics are described by the Landau-Lifshitz equation 
for the constant-length magnetization vector, have been 
studied in details. From a microscopic point of view de- 
scription of such magnets is based on the Landau-Lifshitz 
equation for a unit (normalized) magnetization vector m, 
m 2 = 1, see Refs. [l. 2. 4. Basically, for antiferromagnets 
one can use a set of two equations for magnetizations of 
sublattices, which are unit vectors mi and rri2, or, that 
is more convenient, their irreducible combinations 

m = (mi + m 2 )/2, 1 = (mi - m 2 )/2, (1) 

which are bound by constraint 

(m, 1) = 0, m 2 + l 2 = 1. (2) 

These variables naturally reflect the symmetry inherent 
to antiferromagnets, regarding sublattices rearrangement 
and they are convenient for presentation of phenomeno- 
logical energy of antiferromagnet. However, the grow- 
ing of the number of variables essentially complicates the 
analysis, and within the framework of this a ppro ach a 
few works have been done, we point out Refs. IMSO! 



A considerable progress in study of non-linear dynam- 
ics of antiferromagnets has been reached after obtain- 
ing of so-called a— model, which presents a dynamical 
equation for the antiferromagnet vector 1 see for review 

Refs. Himmmiiiii!. While dedudn s the m ° dei ° ne 

considers that m is small, m 2 <C 1, where m is a slave 
variable and is determined by the vector 1 and its time 
derivative dl/dt. er— model equations can be derived ei- 
ther directly from the Landau-Lifshitz equations for sub- 
lattices magnetizations, 16 ' 17 or phenomenologically, by 
account taken of symmetry considerations. 18 It is a com- 
mon belief that description of nonlinear dynamics of an- 
tiferromagnets within a— model has the same level of uni- 
versality as within the Landau-Lifshitz equations for fer- 
romagnets. At least it is considered to be true for low- 
frequent dynamics in the longwave approximation. 

It is worth noting, the transition to tr-model is not 
connected with any expansion over small amplitudes of 
deviations of the vector 1 from the equilibrium posi- 
tion. Hence, cr— model is highly nonlinear. Since within 
a— model 1 is considered as a unit vector, this model is a 
typical nonlinear chiral model, in which a non-linearity 
is determined a geometric condition 1=1. However, it 
turns out that an isotropic a— model as a nonlinear sys- 
tem is to a certain extent quite "poor" . In particular, for 
a non-localized nonlinear wave of a structure l x + il y = 
lo ■ exp(kr — cot), l z = y/l — Iq — const, where the wave 
amplitude Iq < 1 can be not small, the frequency lj for 
a given "wave vector" k is independent on the wave am- 
plitude Iq. As well, in this system there are no traveling- 
wave solitons, which are most indicative nonlinear exci- 
tations. Note that for the case of anisotropic antiferro- 
magnets with a uniaxial or rhombic magnetic anisotropy 
such traveling-wave solitons a re pre sent, they describe 
moving domain walls, see Refs. ll4Ul5l It is interesting to 
sort whether the abovementioned absence of two specific 
nonlinear effects is an intrinsic property of an isotropic 
antiferromagnet or it appeared due to approximations 
done during transition to a— model. 

To answer this question it is necessary to proceed from 
a full system of equations for two vector variables m and 
1, bound by the relation @. Such an analysis, in princi- 
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pie, is considerably complicated as one has to deal with 
four nonlinear equations, rather than two angular vari- 
ables for the unit vector 1 as for a— model. However, we 
can limit our consideration to analysis of some concrete 
class of solutions in order to confirm the presence of soli- 
tons. 

In this article, a class of solutions in a simple model of 
an antiferromagnet with consideration of only isotropic 
exchange interaction is pointed out. In such a solution 
the vector m is parallel to some direction and change 
its length only, while the vector 1 turns around it within 
some plane. It is appropriate to call these solutions as 
"planar" . Within the class of such solutions, consistent 
description of properties of nonlinear waves and soliton 
dynamics is done. 

The article is organized as following. In the Section 
2 a model is formulated and effective equations of spin 
dynamics in terms of m and 1 without application of typ- 
ical for er— model approximations are presented, and the 
integrals motion are obtained. On the basis of these equa- 
tions in the Section 3 the soliton structure is calculated 
for a one-dimensional case. The analysis of the disper- 
sion law of solitons done in the Section 4 demonstrates 
these one-dimensional stable solitons are magnetic analo- 
gies of Lieb states known from one-dimensional Bose gas 
model. 19 Further in the Section 5 two-dimensional soli- 
tons describing magnetic vortices are analyzed. 

II. MODEL, EFFECTIVE EQUATIONS OF SPIN 
DYNAMICS AND CONSERVATION LAWS. 

Dynamical equations for the vectors m and 1 can be 
written as follows 

/ 6W\ f, SW\ 

^H mx ^M lx ^r)< 

„„S1 / SW\ ( 8W\ , n . 

^ = ( lx ^) + ( mx ^r)' (3) 

where W = W[m, 1] = J w{m, \}{d d x/a d ) is the energy 
functional of an antiferromagnet, which is presented here 
for a magnet with a hypercubic lattice with dimension 
d, w — u>{m, 1} is the energy density, which depends on 
the vectors m and 1 and their spatial derivatives. In the 
standard expansion on gradients with account taken of 
Eq. ^ , a general expression for w in the case of a purely 
isotropic antiferromagnet takes the form 

w = JS 2 m 2 + ^A ia 2 S 2 (Vm) 2 + X -A 2 a 2 S 2 (VI) 2 , (4) 

where J is the effective homogeneous exchange constant, 
the parameters A\ and A 2 are determined by exchange 
integrals within one sublattice and between sublattices, 
respectively, S is an atomic spin and a is the lattice con- 
stant. For this energy, the magnetization vector m equals 
to zero in the ground state. It is worth noting, for the 
model (jlj) with the A\ = the a— model representation 



is exact, while the values of both constants are impor- 
tant for the soliton solutions for antiferromagnet. Below, 
we will not specify relations between the constants A\ 
and A 2 and their connections with some microscopic spin 
model. 

The equations §3§ have the obvious integral of motion, 
the whole system energy E, values of which coincide with 
the value of VF[m, 1] calculated for some concrete solu- 
tion, and the field momentum P, which would be de- 
scribed below. For an isotropic problem the total spin 
value S( tot ) = J Smd d x/a d is also an integral of motion. 
It can be derived from a dynamical equation for the spin 
density m, with usage of energy form ((4]), that gives 

h-^ = -div [S (m x At Vm) + S (1 x A 2 V1)1 . (5) 
at 

This expression determines the conservation law of the 
total spin S(* ot > in differential form. Its analysis allows 
also to point out a concrete exact class of solutions for 
the full set of equations (0). Let the vector m and its 
time and space derivatives at the initial moment of time 
are parallel to some direction, which can be chosen as the 
z axis. The equation ([2]) demonstrate that in this case 
the vector 1 and its derivatives lie in the perpendicular 
(x, y) plane. In virtue of §5§ such geometry remains for 
subsequent moments of time, i.e. dynamical equations for 
an antiferromagnet allow a planar solution in the form of 
m||e z , l_Le 2 . Accounting the constraint © the vectors 
m and 1 can be parameterized by two angular variables, 

m = e z sin fi, 1 = cos fi (e x cos ip + e y sin <p) , (6) 

where e x and e y are unit vectors directed along x and y 
axis, respectively. The initial isotropy of the problem in 
this case manifest itself in arbitrary directions of axis e 2 , 
e x and e y , specific for the planar solution. 

An important characteristic of the planar solution is 
that the system dynamics with new variables/x, ip can be 
described by a simple Lagrangian 

f ddx ( ±* d( P ■ \ 

J -^{~ hS m sln ^- w )' (7) 

where w is the energy density Q presented through an- 
gular variables, 

w = JS 2 sin 2 fi + ^A 2 a 2 S 2 cos 2 fi (Vpf + 

+ ^a 2 S 2 [At cos 2 fi + A 2 sin 2 fi] (V/i) 2 . (8) 

Lagrangian approach allows one to obtain an expression 
for linear momentum of the magnetic excitation P, which 
is a total field momentum of corresponding field, 

P = hS J ^r(V<£) sin^i . (9) 

The dynamical part of the Lagrangian ([7]) and the expres- 
sion for momentum @ contain singularities connected 
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with non differentiability of the azimuthal angle ip. This 
property of the variable ip plays a significant role in de- 
scription of vortices dynamics in ferromagnets.— In our 
case, the presence of this singularity will also manifest it- 
self essentially in description of solitons dynamics, either 
one-dimensional or two-dimensional, see Sections 4,5. 



III. NON-LINEAR WAVES AND 
ONE-DIMENSIONAL SOLITONS. 



Lets us consider dynamics of a simple magnetization 
wave propagating along some direction, say, the x-axis, 
with the velocity v. For such wave, [i = ip = (p(£), 

£ = x — vt. For an analysis of such solutions it is easier to 
start with the spin conservation equations (JSJ) , which can 
be integrated once and then gives an apparent relation 
of ip' — dip/dt; (in this Section, the derivative over £ is 
denoted by prime) and fj, in the following form 



<P 



Hv sin /i. + C\ 
a 2 SA 2 cos 2 \i 



(10) 



where C\ is an arbitrary constant. Using this expression 
it is possible to introduce the Lagrange equation 8L/Sfi 
in the form of the second order ordinary differential equa- 
tion for ji (£). It is easy to demonstrate that this equation 
has the first integral, and for fj, (£) one can obtain a simple 
equation with separating variables. Hence the problem 
allows a general analysis of nonlinear waves depending 
on one parameter, the wave velocity v, and containing, 
in a general case, two arbitrary constants C\ and C 2 . 
The explicit solution of this equation can be presented in 
elliptic functions. 

First of all we are interested in soliton solutions for 
which, far from a soliton, at £ — * ±oo, \x (£) turns zero, 
while tp (£) has constant value. Therefore we consider 
only the case C% = 0. Then the equation for /x(£) in 
soliton solution acquires the following form 



a 2 [Ai cos 2 (i + A 2 sin 2 //] (//) 



/ \ — 



sin /J 2J T7Wa 9~ 



(11) 



Let us discuss properties of such soliton solutions. A 
simple analysis demonstrates that the soliton velocity has 
an upper limit, the value c = 2y 'JAzKS/a, which coin- 
cides with phase velocity of linear excitations (magnons) 
for antiferromagnct. This is a rather natural condition 
for traveling-wave solitons. It is worth noting that c does 
not depend on constant A\, thus, it can be obtained in 
the framework of a— model. However soliton states exist 
only at A\ > 0. The latter is a formal confirmation of 
the fact that for their analysis one should go beyond this 
model. 

The soliton solution of this equation can only be writ- 
ten through elliptic functions. The structure of the pla- 
nar solitons in antiferromagnets, as well as the energy 



dependence on the soliton velocity, is quite common to 
that for solitons in spin nematic stated Hence, we will 
not discuss it in details and limit ourselves with its qual- 
itative analysis. First of all, the form of the solution 
depends highly on the value of soliton velocity v. If the 
velocity v is nearly c, the soliton amplitude /^ ma x is small, 
proportional to \Jc — v. The maximal value of ^t max is 
reached at the zero soliton velocity. 

As follows from the equations (fl"0| and (fTTj) , the values 
of p at the right and left of the soliton differ by a certain 
value Aip. In the case A\ = A 2 the value of Atp = ir and 
it is independent of the soliton velocity. For any other 
relation between Ai and A 2l this limit value A(p = tt 
appears at zero soliton velocity, but Atp < n for v =/= 
and it vanishes at v — > c. In principle, all these features 
are common to that for a so-called rotary waves for easy 
plane ferromagnets, see for review Refs. [HH, or the so- 
called dark solitons, which are well known in nonlinear 
optics^ 2 - 

The energy of a soliton is one of most important soliton 
characteristics. Using Eqs. fIU|) and (fTTj) . the energy 
density w {5} can be easily present through the function 
only. It is convenient to write down the soliton 
energy E as a definite integral over fj, from fi = till 
the maximal value /Lt max - Again, the explicit value of 
this integral can be written through a simple but long 
combination of elliptic integrals only. The exception is 
the limit case A\ = A 2l for which the explicit form for 
soliton energy as a function of its velocity can be written 
as a simple square root dependence, 



E = E \ 1 



(12) 



where c is the spin wave speed, E — 2aS 2 \/2 J A is the 
maximal soliton energy, corresponding to the zero soliton 
velocity v — 0, in the case A~A\ = A 2 . 



IV. SEMICLASSICAL QUANTIZATION OF 
ONE-DIMENSIONAL SOLITONS. 



The soliton energy E and momentum P are the for 
are the most natural soliton characteristics and the de- 
pendence E(P) is the basis for their semiclassical quan- 
tization.— i 2 - Within this approach, E(P) dependence can 
be considered as a dispersion n law for quantum nonlin- 
ear elementary excitations that are described by solitons. 
Usually, this dependence, which is found from classical 
solutions, well reflects the properties of the correspond- 
ing quantum results. 

As has been noted above, the energy is maximal for a 
stationary soliton with v — 0, and it vanishes at v — -> c. 
The concrete dependence can be easily found by numeri- 
cal estimates of corresponding integral, see Ref. [2~D. Con- 
cerning soliton momentum, the situation is not so easy. 
It is worth noting, the equation (TIT))) gives dp/d^ = 
at v — and C\ =0, that formally means zero value 
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of momentum. On the other hand, for any v ^ the 
soliton momentum P(v) is finite, and the limit value 
of the function P(v) at v — > is also finite. For ex- 
ample, for simplest case A\ = A2 one can easily find 
P = (hS/a) ■ arccos (v/c), that gives P — ► ±irhS/2a at 
v — > ±0. Combining this dependence with Eq. (|12|) . one 
can present the dispersion relation for this particular case 
as a periodic function, 



E = E 



nP 
2~Fb 



, Po 



2nKS 



(13) 



with universal period Po- The question appears, whether 
or not these features, the periodicity of the dispersion 
relation and the value of period are model independent. 

In principle, this problem can be overcame by de- 
tail investigation of the behavior of the soliton solution 
at small velocities, see Refs. OH f° r more details. On 
the other hand, it is useful to present a general model- 
free discussion, as it has been done for domain walls in 
ferromagnets i 23 ' 24 Let discuss this problem in more de- 
tails; moreover, it will be useful for the description of 
dynamical properties of vortex-like two-dimensional soli- 
tons. 

Indeed, according to Eq. ([9]), the soliton momentum 
contains a singularity related to the presence of the gra- 
dient of the azimuthal angle (p. Such singularity is an 
internal property of the Lagrangian, see Eq. ([7]) . It be- 
comes clear if we parameterized the spin variables of the 
planar solution through a three-dimensional vector R, 
R = (X, Y, Z) — (m, l x , l y ), whose components repre- 
sent nontrivial variables for the planar solution, namely, 
a magnetization m = m z and two non-zero projections 
of the vector 1. Then the density of the dynamical part 
of the Lagrangian ([7]) can be written as 

where the vector A has a singularity along the Z— axis. 
This Lagrangian coincides with that for a charged par- 
ticle with the coordinate R in a magnetic field with the 
vector potential A. This representation also holds true 
for a ferromagnet in terms of the Landau-Lifshitz equa- 
tion; however, expressions for A in these two cases are 
different. We can readily show that, although the ex- 
pressions for A = A(R) arc different for the cases of an 
antifcrromagnetic planar solution and a ferromagnet, for 
Eq. (TTJ) we have B = rot A = HSR/aR 3 Thus, as in 
the case of a ferromagnet, Eq. (|14[) describes the vector 
potential of a magnetic monopole located at the origin. 
Therefore, the expressions for a momentum P of one- 
dimensional soliton can be obtained by the substitution 
dR/dt — > — vdRd^; it can be reduced to the same form 
as for a soliton in a ferromagnet by gauge transformation. 
We then can use the same method as in Refs. I23ll2~4l 

The formula for the one-dimensional soliton momen- 
tum P = J A(R)dR, contains a singularity and is not 
invariant with respect to the gauge transformations of 



the vector potential A. However, it is important that 
the vector B does not contain singularities on the sphere 
R 2 = 1. Whence, it follows that the difference in the mo- 
menta of two different soliton states is a gauge-invariant 
quantity. Indeed, every soliton (e.g., solitons with differ- 
ent velocities) can be associated with a trajectory con- 
necting certain points R( _ ) and R( + ' lying in the equator 
of the sphere R 2 = 1 (circle Z = or m = 0). In this 
case, the momentum of this soliton is specified by the 
integral J Ac?R over this trajectory going from the point 
R(~) to the point R( + ). Although different solitons (e.g., 
solitons with different velocities) have different values of 
the variable tp gat infinity, all of them have to = at 
infinity; that is, they finish at the equator of the sphere 
R 2 = 1. In this line, the integrand is exactly zero; there- 
fore, the ends of the illustrating trajectories of two soli- 
tons that finish at different points in the great circle can 
be connected by a segment lying in this circle and can 
be considered to be closed. It is clear that the differ- 
ence in the momenta of the two solitons is determined 
by the integral over the closed contour § AdR bound by 
the trajectories describing these solitons. According to 
the Stokes theorem, this integral can be written as a flux 
of a vector B = rotA through the surface enclosed by 
this contour. Therefore, the difference in the momenta 
of two soliton states AP can be represented in the gauge 
invariant form 



AP 



HS 
a 



BdS= 



ns 



cos n dfi dip (15) 



Here the variables tt/2 — /1 and ip can be considered as the 
standard spherical coordinates for the vector R, and the 
integral is taken over the region on the sphere bound by 
the trajectories corresponding to these two solitons. It is 
natural to choose the equator as the line corresponding to 
P = 0, to which the soliton trajectories tend asymmetri- 
cally as the soliton amplitude decreases; this corresponds 
to E — > and v — > c. The maximum soliton energy cor- 
responds to a trajectory that passes through the "north 
pole" of the sphere; for this pole, we have P = Po/2 
and E = P max - The V(P) and E(P) dependencies are 
then qualitatively restored. Indeed, all trajectories cor- 
responding to a soliton velocity in the range from v = c 
to v — or to a soliton momentum from zero to Po /2 fill 
the gap between these two limit trajectories. Hence, the 
momentum increases continuously when going from the 
trajectory near the equator and when approaching the 
limiting trajectory with ±Po/2. As a soliton trajectory 
moves further in the second half of the upper hemisphere, 
the energy decreases and the momentum increases until 
this trajectory reaches the equator. Here, the energy is 
E = 0, the momentum (with allowance for the choice of 
its reference point) is determined by integral (|15| over 
the entire upper hemisphere, and P = Pq. 

Thus, as for domain walls in a ferromagnet , 23 i 24 a true 
periodic E (P) dependence appears for a planar solitons 
in an antiferromagnets due to the topological properties 
of the Lagrangian. This fact should lead to specific fea- 



5 



tures in forced soliton motion, e.g., to oscillating soliton 
motion under the action of a constant force (Bloch oscil- 
lations) as was discussed in details by Kosevich in Ref.l25l. 



V. TWO-DIMENSIONAL SOLITONS - 
ANTIFERROMAGNETIC VORTICES WITH 
FERROMAGNETIC CORE 



V 




Let us consider the static and dynamic properties of 
two-dimensional topological solitons on the basis of the 
model given by Eq. ((¥]). For two-dimensional planar 
solitons the Lagrange equation for the variable ip takes 
the form 



d/u, 

hS— sin u 
at 



A 2 a 2 V[sin 2 /i(V^)], 



(16) 



In the static case, according to this equation, a two- 
dimensional solution can be taken in the form 



(17) 



where r and x are the polar coordinates in the plane 
of the system and (po is an arbitrary angle. To have 
a continuous distribution of the vectors m and 1, the 
number m should be integer. The structure of the vortex 
core is determined by the function fi(r) for which the 
ordinary differential equation can be obtained 



[l + n 2 sin 2 //] 



sm /i cos fj, 



d 2 [i 1 d[i 
dr 2 r dr 

1 A 2 m 2 



A X T 2 



dfi 
dr 



0, (18) 



k 2 = (A 2 — A\) /Ai, Iq = ay/Ai/2 J is the characteris- 
tic length scale. If the condition A\ = A 2 = A holds, 
Eq. (TT8"]) by substitution /i — > it/ 2 — 9 transforms into 
the equation describing the vortex in easy plane ferro- 
magnet, see Refs. [llH. It is easy to show that at r ^> l 
the quantity \x reaches its equilibrium value \i = 0, and 
the behavior near the coordinate origin is a power law: 
/i(r) — 7r/2 cx r rn . Such power dependence is characteris- 
tic of a out-of-plane vortex in ferromagnets. Thus, at the 
center of the planar antiferromagnetic vortex a nonsin- 
gular saturated core with approximately ferromagnetic 
order is formed, and in the vortex center the magnetiza- 
tion takes its maximal value, see Fig.l. 

It is easy to show that the energy of a planar antiferro- 
magnetic vortex, as well as of other topological defects, 
has a weak (logarithmic) divergence with an increase in 
the system size L, it can be written as 



E 



2 nA 2 S 2 a 2 , / L 
m ■ In — — 



(19) 



where r\ is a numerical factor on the order of unity. Hence, 
the vortex with m = ±1 has the minimal energy, and 
further we will discuss only this case. 



y 



FIG. 1: Schematic distribution of the vector 1 (in-plane arrows 
with wide heads) and the vector m (vertical arrows) in the 
planar antiferromagnetic vortex with the vorticity m = 1. 
The core border, chosen as the line with [i — n/4, is marked 
by the dashed line circle. The outermost circle (formally, the 
circle with r — » oo, with the value of /j, — 0) is schematically 
shown by the dotted line circle. 



It is interesting to compare the energy of this pla- 
nar antiferromagnetic vortex with that for vortices in 
easy-plane antiferromagnets. In principle, planar antifer- 
romagnetic vortices contains a ferromagnetic core with 
almost parallel sublattice magnetizations mi and rri2. 
On the first glance, this costs too much energy com- 
paring with that for easy-plane antiferromagnetic vor- 
tices. But this energy difference enters the logarithmic 
multiplier, see Eq. (fT9"|) . Thus, this difference is unim- 
portant for many physical applications; for example, the 
only logarithmic dependence of the energy on the system 
size is manifesting the temperature of the Berezinskii- 
Kosterletz-Thouless transition in two-dimensional sys- 
tems. Thus, both kinds of vortices can be important for a 
description of such transitions for real antiferromagnets. 

Let us describe dynamic properties of the planar an- 
tiferromagnetic vortex, which are also nontrivial. In the 
framework of the a— model, the solution describing any 
soliton freely moving with a velocity of v < c can be ob- 
tained from the known immobile solution by the Lorentz 
transformation with the chosen speed c. However, the 
a— model is inapplicable for the planar antiferromagnetic 
vortex considered above. Analysis shows that the motion 
of the planar antiferromagnetic vortex is possible only 
against the background of "spin flux," i.e., a nonzero 
value of Viy9 — k at infinity. Vortex velocity v and k 
are related as HSv = 2a 2 A 2 ■ k-; this relation can be de- 
rived using the same method as in Ref. [26| for a vortex 
in a ferromagnet. On the other words, far from the core 
of moving vortex the "condensate" is non-uniform, with 
V(f = k cx v = dX/dt. Thus, the total energy of the 
system containing a freely moving planar antiferromag- 
netic vortex diverges as v 2 L 2 , L 2 is a system area, and 
the notion of the local inertial mass losses meaning. This 
property is known for vortices in ferromagnets or super- 
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fluid systems and corresponds to freezing of vortices in 
the condensate, see for review Ref. I3fl5l. 

The problem of the forced motion of the planar anti- 
ferromagnetic vortex can be considered by analyzing the 
field momentum P. Similar to a ferromagnet, Eq. ([8]) 
includes the non-differentiable expression, which leads to 
nontrivial features of the momentum of the topological 
soliton in these systems^ It is most simple to use the 
method proposed in Ref. [27] and to calculate the quantity 
dP/dt in the leading approximation in the vortex veloc- 
ity v. To this end, it is sufficient to use the immobile 
solution given by Eq. (fTU]) with a change of r by f , where 
f = r X(t), X = X(t) = Xe x + Ye y is a coordinate 
of the vortex center. In this approximation, /i = /Lt(f), 
if = mx, r = |r| and \ — arctan[(y — Y)/(x — X)]. Hav- 
ing in mind some general features of the vortex motion for 
the models with gyroscopic dynamics like in Lagrangian 
of Eq. ([7]), let us start with the general form of these 
term as in Eq. (fT4")). not using the concrete form of the 
vector-potential A. 

In the leading approximation on the vortex velocity v, 
the a— th component of the time derivative of the vor- 
tex momentum, dPo /dt with the taken into account the 
conditions dH/dt — —v a (dR/dx a ) can be rewritten as 

dP , a dR z dR 3 fdA, dAA 

— =J d ^di-^ {dR--dR-) (20) 

As for the momentum of one-dimensional soliton, this 
expression contains gauge-invariant quantity B = rotA, 
dAj/dRi — dAi/dRj — Eyfc(rotA)fc, instead of vector- 
potential A as itself. Then the direct calculation yields, 
dPo/dt — G ■ (e z x V). Here the gyroconstant G, as 
well as the linear momentum for one-dimensional solitons 
(|15| , can be presented in the gauge invariant form G — 
fiS J BdS, as a flux of the vector B through the area of 
the sphere R 2 = 1, corresponding to the vortex, that 
gives G = 2nhS/a 2 . 



VI. CONCLUSION. 

Thus, beyond the a— model approximation the 
isotropic antifcrromagnets shows a reach variety of mag- 
netic solitons with non-trivial static and especially dy- 
namic properties. For one-dimensional magnet, soli- 
ton elementary excitations with a periodic dispersion 
law exists. These soliton excitations have common fea- 
tures with the so-called Lieb states,— which are well 
known in many condensed matter models. For two- 
dimensional case, planar antiferromagnetic vortices hav- 
ing non-singular macroscopic core with the saturated 
magnetic moment are found. The dynamic properties of 
these planar antiferromagnetic vortex are also unusual. 
Moving planar antiferromagnetic vortex is subjected to 
the gyroscopic force G ■ [e z , V] , equivalent to the Lorentz 
force for a charged particle in the uniform magnetic field, 
it is well known for vortices in easy-plane ferromagnets 
and superfluid systems, and is observed in experiments 
on the motion of magnetic bubbles and Bloch finest 
In contrast, gyroforce never appears in Lorentz- invariant 
a— model equation; for a usual vortex in an antiferro- 
magnet the gyroscopic force can be induced only by the 
strong external magnetic field and is absent for H= Ch- 
it is worth noting, both these non-trivial dynamical char- 
acteristics, period in dispersion law Pq and gyroconstant 
G, can be written through gauge-invariant expressions of 
the common form. These quantities are independent on 
exchange integrals and depends only on a spin value S 
a single crystal parameter, namely, the interatomic dis- 
tance a. 
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